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We present a semiclassical treatment of a fermionic impurity coupled to Bogolyubov modes of a 
BEC. In the lowest order approximation we find a full solution of an initial value problem, which 
turns out to behave differently in the sub- and supersonic regimes. While in the former case no im- 
purity deceleration is observed, in the latter case non-Markovian dissipation effects kick in resulting 
in slowing down of the fermion. Although this scenario is compatible with the one offered by an 
elementary field theoretical picture at weak coupling, the details of the dynamics turn out to be 
completely different. Fluctuation effects can be taken into account by expansion around the classical 
path. This allows for estimation of the all-important polaron radius as a function of fermion-boson 
coupling strength. As an application we calculate the drag force which is exerted by the BEC on 
the impurity moving with constant velocity v. Contrary to the perturbative result, according to 
which the drag force is ~ u'', it turns out to be proportional to in the semiclassical regime. 

PACS numbers: 03.75.Mn, 71.38.Fp, 67.85.Pq, 78.20.Bh 



INTRODUCTION 

Interacting ultracold gases offer unique opportunities 
for the modeling of realistic quantum systems [T]. Dur- 
ing the recent years it not only became possible to study 
such exotic phenomena as BEC-BCS crossover but also 
modeling of systems, inspired by fractional quantum Hall 
effect might soon be within reach [JHl] • One of very in- 
teresting and old condensed matter physics phenomena, 
which shows up reach physics on the one hand and on the 
other hand might immensely profit from parameter ad- 
justability offered by ultracold gas systems is the polaron 
problem [5j. In its simplest realization it describes elec- 
trons in a semiconductor, which interact with phonons of 
the background lattice. The most important feature of 
the model is a fundamental difference in the properties 
of the system for weak and strong coupling between the 
fermionic and bosonic subsystem. 

It turns out that polaron-like models can very effec- 
tively be simulated by fermionic impurity atoms im- 
mersed into a BEC [SHE]. Alternatively one can couple 
them to a continuum of bosonic excitations of some other 
fermionic species [9] . A direct mapping of most of the re- 
sults from the vintage literature is, however, not possible 
for a number of reasons. First of all, the spectrum of 
bosonic excitations is usually profoundly different from 
those in a semiconductor. The second distinction is the 
coupling mechanism between the subsystems. Some less 
important details such as the presence and nature of cut- 
off parameters might drastically alter the properties of 
a given system as well. But the most important differ- 
ence is at the same time an enormous advantage - it is 
the adjustability of parameters, which is far superior to 
that in solid state realizations. This feature makes it 
possible to advance into new and as yet unexplored pa- 
rameter regions. One such situation is the semiclassical 



limit, in which the impurity can be considered to be an 
almost classical particle. In connection with the polaron 
problem this particular parameter constellation is poorly 
understood and with this Letter we would like to close 
this gap. 

We consider a single fermionic impurity interacting 
with a gas of bosons with masses m (we use polaronic 
units, where the mass of the impurity m/ as well as the 
healing length of the BEC ^ and the reduced Planck con- 
stant h are set to unity, throughout the Letter). The 
Hamiltonian reads 

'P^ \ > 1 ^ ^ 

h kM' .q 

+ ^ViB{q)p{q)al,_^ak', (1) 

h' ,q 

where a^^ and are the creation and annihilation oper- 
ators for a boson with momentum k, p is the momentum 
operator of the impurity, Vb b and Vjb are general boson- 
boson and impurity-boson interaction potentials, and the 
dispersion of the bosons is given by 



/i is the chemical potential. The momentum space den- 
sity of a point particle is related to its position operator 
X via 

p{q) = J d?y6Hv~x)e'ty. (3) 

Since we assume a dilute homogeneous gas, we per- 
form the Bogoliubov approximation by replacing the zero 
mode operators oq, a\ by c-numbers [10]. We further- 
more use contact pseudo-potentials for the interactions 
VBB{q) = 9bb, ViB{q) = giB-, where the constants qbb 
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and giB are related to the s-wave scattering lengths of 
the true potentials via the Lippman-Schwinger equation. 
We then obtain a variant of the Frohlich Hamiltonian 

mm 

k k 

For reasons that will become clear below, we restrict the 
dispersion relation of the Bogoliubov quasiparticles 
and the interaction structure Vk to the expressions valid 
for momenta that are smaller than the healing length in 
the BEC. Therefore, we have 

efc = cfc, Vk = Vak, (5) 

where c = l/{y/2ni) is the sound velocity in the BEC and 
the dimensionless coupling constant a can be related to 



the density n of the bosonic gas and the impurity-boson 
interaction constant via a = ng]g/\/2. 

Derivation of the classical equations of motion (EOMs) 
can be performed using the standard procedure after in- 
tegrating out the environmental (BEC) degrees of free- 
dom, see e. g. pjj. However, in order to preserve the 
convergence of the integrals within the employed approx- 
imation it is necessary to employ a cutoff scheme. For 
reasons which become clear below we use a Gaussian 
^ g-fc /'=c[5n]. The resulting EOM in three dimensions 
is given by 



Xi{t) = - I dw^ijit ~ u,r)xj{u) (6) 
Jo 

with a generalized damping kernel 



(e-'^c'- + l) (c^s^fc^r- + ky) - (e'<- _ l) cs (l + ky) 
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(e^^'^''- - 1) (3cs (4 + 2ky + ky) - kych^) - (^e^'''> + 1) {3kys\ (2 + ky) + ky) 



(7) 



where r = x{t) — x{u). The dependence of the rhs of 
Eq. ^ on r is a major structural difference to the case 
of linear coupling to the environment as in the Caldeira- 
Leggett model. This results in a dissipation dynamics 
which is completely different. The influence of such non- 
Markovian dissipation scenario is weak for short time 
scales and much of the dynamics can be described by 
an effective Caldeira-Legget model. However, the devia- 
tions are drastic for intermediate and long times. In this 
case the solution of the initial value problem can only 
be found numerically [T3]. We solve Eq. (|6| for the ini- 
tial value x(0) = and i(0) = vq. As long as vq < c 
(subsonic initial velocity) the difference to a solution in 
which r = is set on the rhs of Eq. ([7| throughout is min- 
imal, see Fig. [l] On the contrary, in the supersonic case 
Vq > c the non-Markovian dissipation leads to drastic 
discrepancies. In the transient regime we observe dis- 
tinct damped oscillations which can be connected with 
the appealing physical picture emerging in the elemen- 
tary self-consistent harmonic approximation bj, accord- 
ing to which the polaron cloud is modeled as an addi- 
tional mass on a spring attached to the impurity. An im- 
portant point is the fact that both approximations cover 
the nonperturbative regime. 

Despite its efficiency the numerical procedure is not 
able to produce reliable results for the long-time asymp- 



totics of the particle. In order to obtain it one has to rely 
on simplifications of the basic equation (|6]). As we have 
shown above, the memory-like dependence of the rhs on 
the details of the particle trajectory at previous times via 
r is an obstacle. A major simplification is achieved by re- 
alization that r can influence the dynamics only for small 
(t — u), thus we can expand r(t,u) — v{t){t — u) + . . . 
and require v(t) to be self-consistently calculated. Then 
everything simplifies considerably and for the long-time 
asymptotics one obtains the following EOM for the com- 
ponent of the displacement along the initial velocity Vq: 



xit)-\c-x^{t)\) 



{c^x{t))x{t)^ 



(8) 



The rhs is zero for v{t) < c and cx v{t) ^ otherwise. 
This means that the solution to the initial value problem 
x(0) = 0, i(0) = vo with Eq. (fsl) is 



xit) 



Vot 



4/3 



if x[t) < c 
otherwise 



(9) 

Especially the deceleration process of an initially su- 
personic impurity shows up an interesting feature. As 
long as v{t) > c the particle slows down and just above 
the threshold c the velocity decrease follows a power law 
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FIG. 1: (Color online) Results of the numerical solution of 
Eq. Q at short to intermediate timescales. The lower curves 
correspond to an initial condition vq = 0.5c whereas the upper 
curves represent the supersonic case vq = 3c. Dashed lines are 
computed in the assumption r = 0, see text. 
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FIG. 2; (Color online) Asymptotic solution to the initial value 
problem with vo — 2c. The main graph shows the evolution 
of the coordinate. The velocity (inset) decays according to 
a power law with exponent 1/3, until it reaches a value of c, 
after which it remains constant and the particle moves freely. 
The time variable is gauged in such a way that the transition 
from the supersonic to the subsonic regime occurs at t = 0. 



with the exponent 1/3, see Fig. [2] [^T After v drops be- 
low c it does not change any more. From the qualitative 
point of view this behaviour is compatible with that pre- 
dicted by the perturbative calculation in a. However, the 
change of velocity is smooth in the latter case as opposed 
to our result, in which the time derivative of v shows a 
jump. 

As we have already mentioned above can be con- 
nected to the conventional polaron radius. To see this, we 



note that when we take into account quadratic fluctua- 
tions around the classical path x{t) (see e.g. [Ml [15]), 
Eq. ^ becomes a generalized Langevin equation in 
which a noise term £^i{t) has to be added to the rhs. ^{t) 
is a Gaussian stochastic process with average {£,i{t)) = 
and state-dependent autocorrelation function 



(10) 



where Gk{t — u, /?) is the real-time propagator for a boson 
of momentum k at inverse temperature /3. As required by 
the fluctuation-dissipation theorem, in the limit /3 — > cx), 
the rhs reduces to jij{t — u, r)/2. When taking averages 
over realizations of the noise of the full Langevin equa- 
tion, ^i{t) no longer explicitly occurs on the rhs of the 
EOM, but we encounter averages of the form 



(11) 



where we assumed that for every t, u, the random vari- 
able r = x{t) — x{u) can be appropriately approximated 
by having Gaussian statistics. Therefore, from Eq. (Ill, 



we see that a natural cutoff for the fe-sums is provided by 
the mean square of r. In a stationary situation, we would 
expect that (r^) - (r)^ ^ {x{t)x{u)) - {x{t)){x{u)), 
which for u ^ t can be considered to be a measure for 
the squared polaron radius Tpd. Therefore, in this ap- 
proximation 



6/r2 



pol ' 



(12) 



The polaron radius can now be calculated self- 
consistently from a functional expansion of the Langevin 
equation in powers of the noise correlation function. In 
agreement with earlier results based on different meth- 
ods [7], for a coupling constant a of the order of unity, 
we find a polaron radius of the order of the healing length, 
justifying our use of the low-momentum approximations 

We believe, that the phenomena described above can 
be very conveniently measured in state-of-the-art experi- 
ments on ultracold fermion-boson mixtures. We envisage 
a dilute fermionic system being trapped in an optical lat- 
tice, which, in turn, is immersed into a BEG as is done in 
e. g. Refs. [Tni[IZ]. Assuming the interwell fermion hop- 
ping to be completely suppressed (for instance by making 
the lattice very deep) the lattice is then moved through 
the BEG at constant velocity. The fermions are then 
subject to a drag force, which can be measured as a pop- 
ulation imbalance of the low-lying fermion states in the 
lattice by means of Brillouine zone mapping technique 
[18) . Alternatively, after the transients die out one can 
switch off the trapping potentials and infer the force from 
the velocity distributions of freed particles. 

From the theoretical point of view such an experiment 
is very convenient since the particle kinematics is fully 
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FIG. 3: (Color online) Asymptotic behavior of the drag force 
for t — )■ oo as a function of the drag velocity vq. In the 
subsonic regime, it is zero as expected. It jumps to finite value 
just above vo — c, subsequently decaying like 1/vq. Inset: 
drag force as a function of time for drag velocites vo — 0.7c 
(solid line), vo = c (dashed line), vq — 1.2c (dashed-dotted 
line) and vq — 2c (dotted line). 



known. Assuming the impurity/polaron to be a classical 
particle the drag force is precisely given by the rhs of 
Eq. (|6| under the condition r{t,u) = vo(i — u), where 
Vq is the drag velocity. The result of calculations is 
shown in the inset of Fig. ([s]) . While in the subsonic case 
the asymptotic value of the force is zero, in the super- 
sonic case some finite value is reached in the stationary 
situation [22]. Very interesting is the behavior of the drag 
force at t — > c» as a function of vq , which is given by the 
rhs of Eq. see Fig.jsj It follows a 1/uq law just beyond 
the sound velocity with a finite value just above it. This 
is very different from the perturbative result ~ Vq of the 
field theoretical computation, see e.g. [12]. It would be 
very interesting to observe these effects in experiments. 
One possible realization could be ^Li fermionic impuri- 
ties immersed into a BEC of ^■^Na [7]. In this setup the 
typical values for the force evaluate to Fao ~ 10~^® N in 
conventional units. This corresponds to an impurity ac- 
celeration of the order of 1 m/sec^ and therefore seems to 
be well within reach of experimental observation [Tj [17] . 

To conclude, by means of the semiclassical approxima- 
tion and its extensions we have investigated the dynamics 
of an impurity in a BEC and interacting with its Bo- 
golyubov modes. It turns out to crucially depend on the 
velocity of the wave packet and is completely different 
in sub- and supersonic regimes. The resulting mechan- 
ical forces can be estimated and we expect them to be 
measurable in the upcoming experiments. 

The authors would like to thank S. Maier, 
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A classical point like particle does not need a cutoff, of 
course, and summations over all k are necessary. The 
proposed procedure is nothing but a mathematical for- 
mulation of the assumption that the impurity has a finite 
dimension ~ fc^^. 

This power law is different for systems of different dimen- 
sionalities. 

Due to a factorized system preparation at t = and 
resulting polaron cloud formation process the force at 
intermediate times is finite. It roughly corresponds to the 
situation in which the impurity is 'shot' into the BEC 
from the outside. 



